PRO DIAND3P
C=[0,1,-1,-1,1,-1,1,1,-1,1,1,1,-1,1,-1,$
   1,1,1,-1,-1,-1,-1,-1,-1,1,1,-1,-1,1,$
   1,-1,-1,-1,-1,-1,-1,1,1,1,-1,1,-1,1,$
   1,1,-1,1,1,-1,1,-1,-1,1]
A=INTARR(53,53)
A=(C#C+1)/2
A[*,0]=0
A[0,1:52]=1

M=INTARR(95,95)
M[0:47,0:47]=ROTATE(A[0:47,0:47],2)
M[48:94,0:47]=ROTATE(A[1:47,0:47],7)
M[0:47,48:94]=ROTATE(A[0:47,1:47],5)
M[48:94,48:94]=A[1:47,1:47]


S=INTARR(189,189)
S[25,27]=4
S[65,70]=8
S[50,63]=3

B=RANDOMU(1,189,189,POISSON=1)
D=CONVOL(S,M,CENTER=0,/EDGE_ZERO)+B


WINDOW,0,Xsize=189,Ysize=189
TVSCL,D
ZOOM,FACT=3,XSIZE=567,YSIZE=567



END
